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ABSTRACT 

In an effort to better understand the evolution of composite supernova remnants (SNRs) and the 
eventual fate of relativistic particles injected by their pulsars, we present a multifaceted investigation 
of the interaction between a pulsar wind nebula (PWN) and its host SNR G327. 1-1.1. Our 350 ks 
Chandra X-ray observations of SNR G327. 1-1.1 reveal a highly complex morphology; a cometary 
structure resembling a bow shock, prong-like features extending into large arcs in the SNR interior, and 
thermal emission from the SNR shell. Spectral analysis of the non-thermal emission offers clues about 
the origin of the PWN structures, while enhanced abundances in the PWN region provide evidence 
for mixing of supernova ejecta with PWN material. The overall morphology and spectral properties of 
the SNR suggest that the PWN has undergone an asymmetric interaction with the SNR reverse shock 
(RS) that can occur as a result of a density gradient in the ambient medium and/or a moving pulsar 
that displaces the PWN from the center of the remnant. We present hydrodynamical simulations 
of G327. 1-1.1 that show that its morphology and evolution can be described by a ~ 17,000 yr old 
composite SNR that expanded into a density gradient with an orientation perpendicular to the pulsar’s 
motion. We also show that the RS/PWN interaction scenario can reproduce the broadband spectrum 
of the PWN from radio to 7 -ray wavelengths. The analysis and modeling presented in this work have 
important implications for our general understanding of the structure and evolution of composite 
SNRs. 


1. INTRODUCTION 

Composite supernova remnants (SNRs) are a class of 
SNRs for which we observe both a synchrotron-emitting 
pulsar wind nebula (PWN) generated by the wind of a 
highly magnetized, rapidly rotating pulsar, and a shell 
from the supernova (SN) blast wave expanding into the 
ambient interstellar medium (ISM). The PWN is char- 
acterized by an inner shock where the wind flow is ter- 
minated as it enters the expanding nebula, and an outer 
shock at the SN ejecta interface. Its spectrum is deter- 
mined by the evolving population of injected particles 
that spiral in the magnetic field and emit synchrotron 
radiation. This rotating magnetic field also results in jet 
and torus structures that have now been observed in a 
number of systems (review by Gaensler & Slane 2006; 
Kargaltsev & Pavlov 2008) . The interaction of the PWN 
with its host SNR significantly modifies the relativistic 
particle population that eventually escapes into the ISM, 
and it may also give rise to the low-energy particle excess 
that produces 7 -ray emission through inverse-Compton 
scattering (e.g., Abdo et al. 2010; Slane et al. 2010, 2012; 
Temim et al. 2013). This interaction also reveals infor- 
mation about the evolution of the PWN and SNR, the 
properties of the SN ejecta, progenitor star, and the am- 
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bient medium into which the SNR expands. For these 
reasons, the study of composite SNRs, particularly in the 
late stages of evolution when the SN reverse shock (RS) 
interacts with the PWN, is of considerable importance. 

G327. 1-1.1 is a composite SNR that appears to be in 
the late stages of evolution, with a PWN that has likely 
undergone an interaction with the RS (Temim et al. 
2009). At these late stages of composite SNRs’ evolu- 
tion, the RS propagates back into the cold SN ejecta 
in the SNR interior and reaches the boundary of the 
expanding PWN. Since the density of the medium sur- 
rounding a core-collapse SN is unlikely to be uniform, 
the RS typically returns preferentially from the direction 
of higher ambient density. This, along with any offset 
of the PWN due to the motion associated with the pul- 
sar itself, results in an asymmetric PWN/RS interaction 
that disrupts the nebula and gives rise to complex struc- 
tures and mixing of the SN ejecta with PWN material 
(Blondin et al. 2001; van der Swaluw et al. 2004). The 
low-energy particles, whose synchrotron losses during the 
crushing phase have not been extreme, persist as a relic 
radio nebula, while fresh energetic particles injected from 
the pulsar begin to reform an X-ray PWN that is dis- 
placed from the radio emission due the passage of the 
reverse shock and/or the pulsar’s motion. Such mor- 
phology has been observed in composite SNRs, such as 
G327.1-1.1 and MSH 15-56 (Temim et al. 2009, 2013). 
In cases where the pulsar’s motion through the SNR be- 
comes supersonic, the X-ray nebula deforms into a bow 
shock with a cometary morphology, as observed in the 
Mouse (Gaensler et al. 2004). 

SNR G327. 1-1.1 contains a bright central PWN whose 
structure is complex in both radio and X-ray bands. It 
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Fig. 1. — A false color Chandra X-ray image of G327.1-1.1 with the 0.75-1.45 keV emission in red, 1.45-2.58 keV emission in green, and 
2.58-7.0 keV emission in blue a . The Molonglo Observatory Synthesis Telescope (MOST) 843 MHz radio contours (Whiteoak & Green 1996) 
are overlaid on the left panel, with the inner contours outlining the radio relic PWN, and the outer counter the SNR shell, approximately 
17'in diameter. The right panel shows the regions used for the spectral extraction; the white dotted regions correspond to those listed in 
Table 1 and shown in Figure 3, while the numbered regions correspond to those listed in Table 2 and shown in Figures 4 and 5. 
a http://chandra.harvard.edu/photo/2014/mshllg327/ (NASA/CXC/GSFC/T.Temim et al.) 


shows evidence for all of the properties described above; 
an SNR shell, a relic radio nebula resulting from an asym- 
metric RS interaction, a newly reformed X-ray PWN that 
appears to be deforming into a bow shock due to the pul- 
sar’s rapid motion (Temim et al. 2009), and TeV 7-ray 
emission that spatially coincides with the PWN (Acero 
et al. 2012). There is still no evidence for pulsed emission 
from the putative pulsar that is forming the PWN. The 
previous 85 ks XMM-Newton and 50 ks Chandra obser- 
vations revealed the PWN’s unusual structures, such as 
the prong-like features extending into what appeared to 
be a bubble in the SNR. interior. While these observa- 
tions led to the determination of the SNR’s overall prop- 
erties and basic structure, many questions still remain 
about the nature of this system and its evolution. We 
still don’t understand what causes the PWN’s unusual 
morphology, what is its current evolutionary stage, and 
how the injected particle population is altered by the 
reverse shock interaction. In this paper, we present a 
study of the deep 350 ks Chandra observations of the 
composite SNR G327. 1-1.1, with the goal of understand- 
ing its physical properties and evolutionary stage. We 
also present a 2-dimensional (2D) hydrodynamical (HD) 
model of G327. 1-1.1 that provides insight into the evo- 
lution and origin of the morphology of this remnant and 
composite SNRs in general. 

The paper is organized as follows. Sections 2 and 3 
describe the observations and basic analysis of the data. 
Section 4 discusses the general morphology of the SNR 
and the PWN structures revealed by the deep Chandra 
observation. Secion 5 describes the spectral fitting re- 
sults and discusses the spectral properties of the SNR 
shell and the PWN. Section 6 presents the HD model for 
the evolution of G327. 1-1.1, simulated as an SNR blast- 
wave expanding in an ISM density gradient and contain- 
ing a PWN produced by a rapidly moving pulsar. Sec- 


tion 7 presents our modeling of the broadband emission 
of G327. 1-1.1 from radio to 7-ray wavelengths. Section 8 
summarizes the results on the evolutionary state of the 
SNR and some of the outstanding questions. The con- 
clusions are summarized in Section 9. 

2. OBSERVATIONS AND DATA REDUCTION 

A deep observation of the SNR G327. 1-1.1 was car- 
ried out with the Advanced CCD Imaging Spectrometer, 
ACIS-I, on board the Chandra X-ray Observatory for a 
total exposure time of 350 ks. The observations were 
taken in the VFAINT mode, on 2012 May 22, 24 and 
27, under the observation IDs 13767, 13768, and 14430, 
and correspond to exposure times of 143, 170, and 37 ks, 
respectively. The standard data reduction and cleaning 
were performed using the chandra-repro script in CIAO 
version 4.5, resulting in a total clean exposure time of 
337.5 ks. 

The merged, exposure corrected and adaptively 
smoothed false color image of G327. 1-1.1 is shown in Fig- 
ure 1, with the 0.75-1.45 keV shown in red, 1.45-2.58 keV 
in green, and 2.58-7.0 keV emission in blue. The hard 
X-ray emission shown in blue originates from the PWN 
material, while the soft X-ray emission in red shows the 
spatial distribution of the thermal emission. The overlaid 
radio contours from the Molonglo Observatory Synthe- 
sis Telescope (MOST) 843 MHz observation (Whiteoak 
&; Green 1996) show the outline of the SNR shell (17' 
in diameter) and the inner relic PWN. The right panel 
of Figure 1 shows the regions used for the spectral ex- 
traction. Spectra were extracted and response files gener- 
ated for each observation ID individually, using the CIAO 
specextract script, producing a set of three spectra and 
associated files for each spectral region. Since the obser- 
vations were carried out close in time and have the same 
roll angle, the spectra and response files were merged to 
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Fig. 2. — Top panel shows the unsmoothed Chandra ACIS image 
of the PWN region (2'6 across), corresponding to region 11 in 
Figure 1 (right). The black box is 20" in size and indicates the 
zoomed-in region of the compact source shown in the bottom panel. 
The bottom panel’s left image shows the unbinned image of the 
compact source, the middle panel shows the residual image after 
the subtraction of a point source at the location of the black “x” , 
and the right panel shows the same residual deconvolved by the 
task arestore. The scale bar is the units of counts. 

produce a single spectrum for each extraction region. 


3. ANALYSIS 

In order to study the properties of the various 
structures observed in the PWN and the spatial 
variations in the thermal emission across the SNR 
shell, we fitted spectra from 23 different source 
regions, shown in Figure 1. Since the SNR shell 
covers most of the field of view, we used the ACIS 
blank-sky background files adopted to our observation 
(http: / / cxc.harvard.edu / ciao / threads / acisbackground/) 
to extract background spectra from each extraction 
region in Figure 1. We used the high energy data 
from 10-12 keV to compare the particle background in 
each ACIS-I chip in our observations to the particle 
background in the blank-sky data. We found that 
the particle background in the blank-sky data was 
~ 18% higher, and we therefore scaled the exposure 
time of the blank-sky data accordingly. We extracted 
backgrounds from two small regions in our observation 
of G327. 1-1.1 that were free of SNR emission (regions 
marked “24” in Figure 1), in order to compare them 
to the corresponding scaled blank-sky background 
spectra. We found that the scaled blank-sky spectra 
described the background very well, except for a residual 
thermal excess between 1-2 keV in the background of 
our observation. Since the blank- sky data only account 


for background emission for high Galactic latitudes, the 
presence of an additional thermal excess in the vicinity 
of G327. 1-1.1, which is located in the Galactic plane, is 
not surprising. We modeled this excess simultaneously 
in the two background regions in Figure 1 and found 
that it is well described by an XSVNEI model with an 
absorbing column density Nh — 1 07 x 10 22 cm -2 and a 
temperature of 0.68 keV. We note that these parameters 
are quite different than those for the thermal emission 
associated with the SNR. 

In order to fit the source spectra from each extrac- 
tion region, we first subtracted the corresponding blank- 
sky spectrum, and then included a model component 
for this additional thermal background excess, normal- 
ized by the extraction area for each region. The source 
and background spectra were then fitted simultaneously. 
All regions were fitted by either a power-law model de- 
scribing the synchrotron emission from the PWN, a non- 
equilibrium ionization thermal model (XSVNEI) with so- 
lar abundances Anders Sz Grevesse (1989), or both com- 
ponents for regions where both non-thermal and thermal 
emission was detected. The results of the fits are listed 
in Tables 1 and 2. Table 1 lists the best-fit parame- 
ters for the spectra extracted from the whole PWN and 
SNR shell, shown in Figure 3. These regions are the 
white dotted regions in Figure 1, where the region for 
the SNR shell is the outer circle with the inner PWN 
region excluded. Table 2 lists the best-fit parameters for 
the numbered regions of Figure 1. The spectra and best- 
fit models for all the regions are shown in Figures 4 and 
5. 

4. GENERAL MORPHOLOGY 

The radio morphology of G327. 1-1.1 is that of a classic 
composite SNR; a faint shell ~ 17' in diameter surrounds 
a bright PWN, approximately 5' in diameter, located off- 
center from the shell (Whiteoak & Green 1996). The ra- 
dio contours are overlaid on the Chandra X-ray Image in 
Figure 1 (left). A distinct finger- like structure protrudes 
from the PWN in the northwest direction, possibly sug- 
gesting a picture in which a fast-moving pulsar is moving 
through the SNR, leaving the PWN in its wake. X-ray 
studies have confirmed the composite nature of G327.1- 
1.1, and provided evidence of a compact core that is the 
likely counterpart to the pulsar that powers the observed 
PWN emission (Lamb & Markert 1981; Seward et al. 
1996; Sun et al. 1999; Bocchino & Bandiera 2003; Temim 
et al. 2009). Imaging in the hard X-ray band shows that 
the high energy photons are clearly offset from the main 
radio PWN in the direction of the radio finger, sugges- 
tive of a scenario in which the unusual morphology is 
produced by a combination of the pulsar’s motion and 
the disruption of the PWN by the passage of the RS. 
XMM-Newton observations of G327. 1-1.1 showed a ther- 
mal X-ray shell, coincident with the radio shell emission, 
with brightest emission in the southeast (Temim et al. 
2009). 

The new 350 ks Chandra observation presented in this 
work also shows soft X-ray emission that coincides with 
the radio shell, as shown by the red color in Figure 1. 
This diffuse emission is enhanced in the southeastern 
part of the shell, corresponding to region 18. The X- 
ray PWN is clearly displaced from the geometric center 
of the shell. The deep Chandra images show that the 
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Fig. 3. — Best-fit models to the spectra for the whole PWN and SNR shell regions, the white dotted regions in the right panel of Figure 1, 
where the extraction region for the shell is the annulus that excludes the inner PWN region. For all three panels, the power-law component 
is shown in blue, thermal component in red, and residual background in gray (see Section 3). The left and middle panels show the spectrum 
and best-fit models for the PWN region, which is equally well fitted by a thermal model with enhanced and solar abundances, respectively. 
The right panel shows the spectrum and best-fit model for the SNR shell. The values for the best-fit parameters are listed in Table 1. 


non-thermal emission from the PWN extends beyond the 
previously identified trail that coincides with the radio 
finger (Temim et al. 2009). The full extent of the hard 
X-ray emission from the PWN, shown in blue in Figure 1, 
approximately fills an ellipse with major and minor axis 
length of approximately 7'9 and 4'4, respectively. X- 
ray emission is now seen to completely fill the radio relic 
PWN, outlined in the contours of Figure 1. In addition 
to hard X-ray emission, the inner region of the radio relic 
also contains softer X-ray emission that is likely associ- 
ated with thermal emission. While this emission may 
arise from the swept-up ISM seen in projection, its spa- 
tial correlation with the relic nebula suggests that it may 
originate from the inner SN ejecta that have mixed with 
the PWN material following the passage of the RS, as 
previously observed in systems such as Vela X (LaMassa 
et al. 2008). 

As revealed by the earlier Chandra data, the brightest 
part of the PWN is the compact core at the tip of the ra- 
dio finger which appears to contain the neutron star that 
powers the central PWN (Temim et al. 2009). The core 
is embedded in a bow-shock shaped cometary structure, 
with an extended X-ray trail connecting to the relic neb- 
ula, and prong-like structures that extend into large arcs 
(~ 2. '5 in length) towards the northwest (regions 13 and 
14 in Figure 1) that appeared to form a bubble (Temim 
et al. 2009). The deep Chandra observation provides no 
evidence that the arcs actually connect to form a bub- 
ble, but instead suggests that this “outflow” is an open 
structure. 

4.1. Structure of the Inner PWN 

The 350 ks Chandra observation presented in this work 
provides a much more detailed view into the structure of 
the inner X-ray PWN, contained in region 11 in the right 
panel of Figure 1. The zoomed in image of this inner re- 
gion is shown in Figure 2. A bright compact source, 
~ 10" in diameter, is embedded in a 40" long cometary 
structure whose morphology has been suggestive to re- 
semble a bow shock produced by a pulsar moving super- 
sonically through the SNR interior (Temim et al. 2009). 
Behind the cometary region is a trail of X-ray emission, 
approximately I'A in length, that spatially coincides with 
the radio finger and connects the cometary structure 


TABLE 1 

SPECTRAL FITTING RESULTS 1 


PARAMETER 

PWN 

Shell 

Area (arcsec 2 ) 

207509 

516889 

Counts 

147448 

65420 

N h (10 22 cm- 2 ) 

l-76 +0:11 

o-qq+u:^ 1 

0 - 09 + 9 : i4 

Photon Index 

o-i r+0:05 
z - io -0:04 

-0:02 

0.7Q+0-23 

^■'°- 0:25 

Amplitude (10 -3 ) 

q-ro + 0-31 
O.U°_0;19 

4:46±g;It 

1-4C+0.52 

Fi (10 _11 erg cm -2 s _1 

1:88 

2:24 

0:71 

kT (keV) 

0 . 2 q+0:14 

u - zy -0:09 

0-29 +0:01 

u - zy -0:06 

0‘29 +0:05 

u - zy -0:04 

[Mg]=[Mg]© 

q.Q+2. 5 

05.8_o: 8 

1.0 

1 .q+0.5 

l" 5 — 0:2 

[SiWSiJo 

r.i +5:7 

1.0 

1 .0 + 1.1 

— 0:4 

(10 10 s cm -3 ) 

6:8+S s 

0‘57+ 0:36 

°:33+b : : 14 

Norm. (10 2 ) 

0:83 ±g;g! 

6‘51 + 9:37 

— 0:87 

8:4611!*! 

F 2 (10 _11 erg cm 2 s -1 ) 

4:8 

16 

24 

Reduced 2 

0.9 

1.0 

1.0 


Note. — The best fit parameters for a power-law plus non- 
equilibrium ionization thermal model (XSVNEI). The listed 
uncertainties are 1.6 sigma (90 % confidence) statistical uncer- 
tainties only. The fluxes were calculated in the 0.3-10.0 keV 
range. The amplitude for the power-law component is in the 
units of counts keV - 1 s — 1 (at 1 keV), and the normalization 
of the thermal model is equal to 10 14 n e nH V=4 d 2 cm 5 , 
where V is the volume of the emitting region and d is the dis- 
tance to the SNR. The parameter represents the ionization 
timescale for the best-fit thermal model, and the number of 
counts for each region represents blank-sky-subtracted source 
counts. The fluxes Fi and F 2 are in the units of erg cm -2 s -1 
and they represent unabsorbed fluxes for the power-law and 
thermal components in the 0.3-10 keV band, respectively. 

back to the radio nebula. We refer to this feature as 
the X-ray “trail”. A pair of narrow, well defined prong- 
like structures, ~ 1'8 in length, are seen on each side of 
the cometary structure, parallel with its long axis. These 
prong structures extend into the more diffuse arcs seen 
in Figure 1. 

4.1.1. The Compact Source 

Previous spatial modeling of the compact source itself 
showed that the emission can be described by a point 
source embedded in a more extended, slightly elliptical 
Gaussian component (Temim et al. 2009). When these 
best-fit components were subtracted from the original 
50 ks Chandra image, there was an indication of an asym- 
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Fig. 4. — Best-fit models and residuals for the spectra from regions 1-12 in Figure 1 and listed in Table 2, where the numbers in the 
upper lefthand corners of the panels correspond to the region numbers. For all panels, the power-law component is shown in blue, thermal 
component in red, and residual background after blank-sky subtraction in gray (see Section 3). 
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Fig. 5. — Best-fit models and residuals the spectra from regions 13—24 shown in Figure 1 and listed in Table 2, where the numbers in the 
upper lefthand corners of the panels correspond to the region numbers. For all panels, the power-law component is shown in blue, thermal 
component in red, and residual background after blank-sky subtraction in gray (see Section 3). 
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metric residual that might suggest a more complex struc- 
ture. The new Chandra observation provides ~ 11,000 
counts from the compact source, approximately seven 
times more than in the original observation. Approxi- 
mately 5,000 of these counts are likely from the under- 
lying cometary structure in which the compact source is 
embedded, and ~ 6,000 from the source itself. A spa- 
tial model of a point source embedded in an extended 
component leads to parameters consistent with the orig- 
inal results of Temim et al. (2009). The point source 
contains ~ 580 counts, or ~ 10% of the total compact 
source flux. It is embedded in an extended component 
well described by a 2D Gaussian with a full width half 
maximum (FWHM) of 5'/ 7 ± O'/I and an ellipticity of 

0. 26. ± 0.02, containing ~ 5,400 counts. Even though 
the overall shape of this more extended component can 
be described by the elliptical Gaussian, its emission is 
not entirely smooth and symmetric, but instead shows 
evidence of a more complicated structure. The bottom 
three panels of Figure 2 show the unsmoothed zoomed-in 
image of the compact source on the left, the residual im- 
age of the extended component after the subtraction of a 
point source at the location of the black “x” in the mid- 
dle, and a deconvolved image of this same point-source- 
subtracted extended component on the right. The decon- 
volution was performed with the CIAO Sherpa task are- 
store that uses the ACIS point spread function to recover 
spatial resolution. Both the original and deconvolved 
residual images (middle and right) clearly show struc- 
ture within the compact source that contains the pre- 
sumed pulsar, reminiscent of a torus structure observed 
around many other young neutron stars (Kargaltsev &; 
Pavlov 2008). 

5. SPECTRAL PROPERTIES 

The 50 ks Chandra and 85 ks XMM-Newton observa- 
tions (Temim et al. 2009) showed that the X-ray spec- 
trum of the PWN can be described by a power-law with 
an average photon index of 2. 11 ±0.03 and the shell emis- 
sion by a Raymond-Smith thermal plasma model with 
a temperature of 0.3 keV. There was also evidence for 
steepening of the photon-index with distance from the 
compact source that may be suggestive of synchrotron 
aging. Our deep Chandra observations confirm the global 
spectral properties found in the previous work (see Table 
1), but they allowed us to study the spectral properties of 
the individual PWN structures, and explore the spatial 
variations in the absorbing column density and the ther- 
mal emission from the SNR shell. Based on the best-fit 
parameters for the source spectra extracted from regions 
in Figure 1 and listed in Tables 1 and 2, we make the 
following observations: 

1. The compact source and the cometary structure (re- 
gions 1 and 2) have photon indices that are consis- 
tent within the uncertainties, with T = 1.61 to :07 anc ^ 
T = 1.62lo5o7> respectively. There is no evidence 
for significant steepening in the photon index in the 
cometary structure. 

2. The photon index of the trail steepens from T — 
1.76 ± 0.12 to T = 2.28 ± 0.12, giving a difference 
of AT (5 — 8) = 0.52 ± 0.17, where the region numbers 
(Table 2) for which the difference in the photon index 


was calculated are indicated in the parentheses. As 
will be discussed in Section 6.2, this can be explained 
by synchrotron aging of particles as they travel from 
the pulsar along the extent of the trail, a distance of 
approximately 80" (3.5 pc at a distance of 9 kpc). The 
trail may be interpreted as PWN material left behind 
by the pulsar as it travels through the SNR, along 
with being swept back by the passage of the SN re- 
verse shock that effectively arrives at the pulsar from 
the NW direction, as a result of an east/west ISM 
density gradient (see Section 6). 

3. The prong structures seem to extend southward of the 
compact source, parallel to the sides of the cometary 
structure. More specifically, regions 3 and 9 and re- 
gions 4 and 10 in Figure 1 seem to be part of the 
same structures (also see Figure 2). These struc- 
ture have a photon index of ~1.8. There is no evi- 
dence for spectral steepening along the prongs, with 
Ar(3-9) = 0.12±0.18 and AT(4-10) = -0.05±0.18. 
The prongs extend into more diffuse and extended 
arc-like features (regions 13 and 14), approximately 
155" in length from the northern tip of the prongs. 
Their synchrotron spectrum is significantly steeper, 
with a change in the photon index from the prong to 
the arc of Ar(13 — 9) = 0.58 ± 0.28 in the east and 
Ar(14 — 10) = 0.65 ± 0.23 in the west. The space 
between the arcs (regions 15 and 16) is also filled with 
synchrotron emission with a photon index similar to 
those of the arcs, suggesting that the whole structure 
may be a funnel shaped outflow, with the arcs possi- 
bly being limb brightened edges, rather than separate 
structures. The surface brightness of the PWN be- 
tween the prongs (between regions 9 and 10) appears 
to be lower than in the surrounding nebula, which 
may be hinting that the funnel shaped outflow may 
be hollow. 

4. The steepest synchrotron spectrum in the PWN is ob- 
served in region 12, corresponding to the relic PWN 
observed at radio wavelengths. The best-fit photon in- 
dex in this region is T = 2.6±0.1. Fainter non-thermal 
emission also seems to be present in regions extracted 
from the SNR shell, away from the central PWN (re- 
gions 17, 18, 19, 21, 23). This non-thermal emission 
in the outer regions of the SNR is faint, at the level of 
only ~ 2% of the total PWN emission. It also has a 
photon index that is considerably steeper, leading us 
to conclude that it is most likely a dust-scattering halo 
from the PWN that is expected to have a softer spec- 
trum due to the energy dependence of the scattering 
efficiency (e.g. Bocchino et al. 2005). 

5. The thermal emission both around the PWN and in 
the SNR shell can be characterized by the same aver- 
age temperature of 0.29 keV (see Table 1). Within the 
uncertainties, there is no evidence that there are sig- 
nificant spatial variations in the thermal temperature 
across the SNR. The SNR properties derived using the 
Sedov (1959) blast-wave model and the parameters 
from the best-fit thermal component are consistent 
with those found in Temim et al. (2009); an ambient 
ISM density no = 0.12 cm -3 , a shock velocity v s = 
500 km s -1 , an explosion energy £ 5 i(10 51 erg s x ) = 
0.5, and an age of 17400 yr. The estimates assume an 
SNR distance of 9.0 kpc. 
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TABLE 2 

SPECTRAL FITTING RESULTS 2 



Region 

Area 

Cts 

N h 

Photon 

Amplitude 

kT 


Norm. 

Fi 

F? 

Red. 



(arcsec 2 ) 

(1000) 

(10 22 cm -2 ) 

Index 

(10- 4 ) 

(keV) 

(10 12 s cm -3 ) 

(10- 3 ) 

(10- 

+ 

2 

1 

Compact Source 

84.657 

6.34 

-|-qo+0:08 
1 ‘ y — 0:08 

i. 61 +0:08 

i - D1 — 0:07 

1 05+° : 11 
I.UO— 0:10 




0:45 


0.80 

2 

Cometary PWN 

971.22 

7.75 

1.93 


1-47+0:16 




1:09 



3 

Trail East 

537.42 

2.13 

1.93 


0-44+0:07 

u>44 -0:06 




0:27 



4 

Trail West 

766.56 

3.12 

1.93 


0'61 +009 




0:39 



5 

Trail 1 

424.45 

1.98 

1.93 

1 '76+°' 12 
4 - ‘ 0: 12 

q-qq+0:05 

u.oy_ 0 :05 




0:26 



6 

Trail 2 

588.19 

2.13 

1.93 

1- 95+0- 11 

0'49 + °' 07 

u,4y -0:06 




0:28 



7 

Trail 3 

994.92 

2.99 

1.93 

2‘09 +0:1 ° 
z.uy_ o;io 

fl-78 4 "®'®® 
U - 0:08 




0:42 



8 

Trail 4 

839.48 

2.38 

1.93 

9-98+®- 42 
2 ‘ 28 -0:12 

0-74+0:09 
u - '^-0:09 




0:37 



9 

Prong East 

828.58 

1.66 

1.93 

,•72+0:14 
-*-• ‘ 0: 14 

0‘30 +0:06 

U.OU_ 0:05 




0:27 



10 

Prong West 

971.22 

2.06 

1.93 

V o,-+0:14 

n-44+ 0:08 

u -^*-0:07 




1:09 



11 

Diffuse PWN* 

20007 

27.7 

1.93 

9-, ,+0:04 
— 0:05 

6'91 + 0:37 
Dyl - 0:74 

0‘23 +0:14 

0:05 

0-91 + 0:88 
u - 2li -0:16 

6:Ol! 6 o 

3:68 

17:7 

0.82 

12 

Relic PWN* 

26787 

17.2 

1.93 

2'58+°’ 07 

651+ 0 - 53 

0.23 

0:21 

6:9l! 8 5 

3:14 

20:3 


13 

Arc East 

6147.6 

4.60 

2:08«;i| 

2:30;S;?S 

1 '62+0:52 
1,0 —0:38 



0:81 


0.74 

14 

Arc West 

6226.4 

5.56 

3:Olt2;S 

2:5ot2;;i; 

9-74+0.74 

0-21+®- 33 

0T4 +O ' 22 

04+8 4 

1:3 

102 

0.95 

15 

Inner Arcs 1 

3430.5 

3.67 

2'94 +0 ' 71 

^' y4 -0:68 

2’27+°' 27 

^• z '-0:22 

1 '61 +® :93 
iD1 — 0:45 

0:2912:!! 

0:6912:54 

2:7± 2:1 

0:61 

6.0 

0.88 

16 

Inner Arcs 2 

22014 

7.23 

2-iq+ 0 - 34 

9. C9+0.22 
z - oz -0:20 

2:74+2:22 

0'29 + °' 52 

2:7+ 2:1 

0-1+26 

1:3 

7:2 

0.88 

17 

East Shell 

42053 

4.43 

1'66 +0:46 

1.00—0:28 

9-ac:+0:62 

z.oo_ 0;92 

0:8012;“ 

0'30" 1 ’ 014 

u.^u_ 0:09 

0:08lg;22 

2:2l! 7 7 

0.40 

11.2 

1.13 

18 

SE Shell 

62109 

11.2 

2'00 + °' 17 

z.uu_ 0:08 

9 . 7 o+0:48 

1-47+ 1 - 61 

0-91+ 002 

0:12+2:22 

1 0-0+34.8 

0.72 

57 

1.02 

19 

North Shell 

18470 

3.28 

rt.rti +0: 53 
zzl -0:70 

3‘01 + 0:73 
,3 - ul -0:69 

1 -40+1’ 94 

i ’^ u -0:84 

o:i8lS:25 

0:08lg;2 4 

49122“ 

0.73 

130 

1.00 

20 

NW Shell 

23944 

2.61 

2' 34+°' 18 



0:53+g;!5 

20:5+::: 

0:85+2:51 


1:5 

0.80 

21 

West Shell 

18227 

3.40 

S'OO" 1 ’ 023 

o.uu_ 0 :3o 

1 -afi+®-® 2 
l-o6_ 0^4 

0'34 + ® 22 
u.^_ 0:12 

0-->7+0.05 
U " 5 ' -0:04 

16:41+3 

3:412:2 

0.21 

6:9 

1.12 

22 

SW Shell 

69487 

7.01 

9.7C+O.29 



O.9Q+O.O6 

u - zy -0:05 

l:8+i:o 

ll+ie 


22 

1.01 

23 

NE Shell 

26361 

2.73 

1 '76 +0:23 

1 ' ' D — 0:32 

2‘46 +0:34 

Z.HO_ 0 ;42 

0'50 +<j:3£> 

U.OU— o:i 6 

O-07+0: 19 
-0:08 

15+::: 

0'36+° 52 
u. oo_o:24 

0.24 

0:79 

0.80 

24 

Background 

19611 

1.55 

l'07 + ° 03 
' —0:09 



n-«s+°- 48 
U.68_ 0:05 

49l is 

n-ns+°- 02 

0.08_ 0:02 


0:21 

0.87 


Note. — Best-fit parameters for the spectra extracted from regions shown in Figure 1, using a power- law and/or a non-equilibrium ionization 
thermal model (XSVNEI). Listed uncertainties are 1.6 (90 % confidence) statistical uncertainties from the fit. The Nh value listed for Region 1 

is the result of the simultaneous fit for regions 1-12. The parameters for the thermal components for Regions 11 and 12 were fitted simultaneously, 
as well as the thermal components for Regions 15 and 16. The parameters for the background region 24 is a simultaneous fit to the two background 
regions labeled “24” in Figure 1. The amplitude for the power-law component is in the units of counts keV - 1 s -1 (at 1 keV), and the normalization 
of the thermal model is equal to 10 -14 n e riH V=4 d 2 cm -5 , where V is the volume of the emitting region and d is the distance to the SNR. The 
fluxes Fi and F2 are in the units of erg cm -2 s -1 and they represent unabsorbed fluxes for the power-law and thermal components in the 0.3-10 
keV band, respectively. *The thermal component for the diffuse and relic PWN regions show a [Mg]=[Mg]© overabundance of 2:3ljv. 61 . The 
thermal component in Region 24 is the excess background after blank-sky background subtraction. The number of counts for each region represents 
blank-sky-subtracted source counts. 


The spectrum extracted from the entire PWN region 
is well fitted by both a model with fixed solar abun- 
dances, and a model with enhanced abundances of Mg 
and Si. While the enhanced abundance model pro- 
vides a slightly better fit, the difference in reduced x 2 
values is not significant enough to distinguish between 
the two models. Both spectral fits are shown in Fig- 
ure 3 and parameters listed in Table 1. The enhanced 
abundance model gives best fit values of 3.8+g:| and 
times solar, for Mg and Si respectively, and 
would suggest an interaction of the PWN with SN 
ejecta, or possible mixing of the SN ejecta and PWN 
material following a RS interaction. As mentioned in 
the previous section, the morphology of the soft X- 
ray emission inside the relic PWN, supports the idea 
that the soft thermal emission in this region originates 
from SN ejecta. This is also supported by the fact that 
the best-fit ionization timescale for the solar abun- 
dance model is implausibly low for emission arising 
from swept-up ISM material in an evolved SNR. The 
most likely scenario is that this region contains a mix 
of SN ejecta and ISM/CSM emission that is either pro- 
jected emission from shocked SNR shell, ISM material 


that has mixed into the SNR interior, or both. 

6. The average best-fit column density is Nh = 
l-76io:o7 x 10 22 cm_2 f° r the PW T N region and N H = 
2.22^0:28 x 1022 cm_2 for the SNR she11 ( see Figure 1 
and Table 1). There is an indication of a somewhat 
higher absorbing column density in the western side 
of the SNR shell (regions 14, 21, and 22). This higher 
N h value along the line of sight in the west is not 
necessarily intrinsic to the SNR, but it is consistent 
with a scenario in which the SNR expands into denser 
ambient medium in the west. 

6. HYDRODYNAMICAL MODELING 
In order to understand the morphology of G327.1- 
1.1, we have modified the VH-1 hydrodynamics code 7 
(Blondin et al. 2001) to model the evolution of a com- 
posite SNR that expands into a non-uniform ambient 
density and contains a rapidly moving pulsar. The 2D 
simulation is evolved on a cylindrical grid (r,0), but with 
spherical radial divergence (volume proportional to r 3 ). 

7 http:/ /wonka.physics. ncsu.edu/pub/VH-l/ 
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Fig. 6. — Results of the HD simulation of a composite SNR expanding in an ambient density gradient, and containing a rapidly moving 
pulsar. The figure shows density maps (g cm -3 ) at four different SNR ages. The density is higher in the west and decreases towards 
the east. The pulsar is moving to the north at a velocity of 400 km s -1 . The input parameters for the simulation are consistent with 
observed values for G327. 1-1.1. The resulting morphology at 17000 years closely resembles the morphology of G327. 1-1.1 at radio and 
X-ray wavelengths. 


The model is constructed within the pulsar’s rest frame 
and consists of a stationary, expanding gas bubble with 
an adiabatic index 7 = 4/3 for relativistic gas, contained 
within an SNR modeled by a self-similar driven wave ex- 
panding into an SN ejecta density profile described in 
Blondin et al. (2001). Since the pulsar is stationary in 
the grid’s reference frame, the SNR and the surround- 
ing ISM were given a negative velocity to simulate the 
pulsar’s true motion. This choice of frame is constructed 
to utilize the inner boundary to act as both an energy 
injection point for the pulsar wind and as a hard-surface 
modified to keep large amounts of mass from exiting the 
grid during the collision with the RS. 

The model assumes a time-dependent pulsar spin-down 
power E = £ 0 (1 + «/" r o) -(n+ 1 H n ~ 1) that drives the 
PWN, where tq is the pulsar’s spin-down timescale and 
n is the braking index. The magnetic field and brak- 
ing index are assumed to be constant, as presented in 
(Blondin et al. 2001). The model also includes a gradi- 
ent within the ambient ISM density that can be rotated 
about the grid center to achieve a desired angle between 
the orientation of the gradient and the pulsar’s direction 


of motion. 

6 . 1 . Input and Output Parameters 

The input parameters for the HD simulation are sum- 
marized in Table 3. We used an explosion energy 
of 0.5 x 10 51 erg and an average ambient density of 
0.12 cm -3 , as estimated from the Sedov model and the 
best-fit thermal parameters for G327. 1-1.1 (see Section 
5). We assumed an ISM density gradient that decreases 
from west to east, as suggested by our best-fit measure- 
ments of the absorbing column density Ny\ . The pulsar 
is assumed to be moving to the north with a velocity of 
400 kms -1 (Temim et al. 2009). 

The parameters for which we did not have obser- 
vational constraints were adjusted to produce a mor- 
phology that most closely resembles that of G327. 1-1.1 
at its estimated age of ~ 17000 yr. We assumed an 
ISM density gradient profile described by equation 14 
of Blondin et al. (2001), with x — 1.08, correspond- 
ing to a density contrast of 12.5, and the characteris- 
tic length scale over which the density changes (H) set 
to 5.2 pc. For the PWN, we fixed the braking index 
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TABLE 3 

SNR PROPERTIES & HD MODEL INPUT PARAMETERS 


Parameter 

Description 

Value 

SNR Properties: 

D (kpc)* 

SNR distance 

9.0 

RSN R (pc)* 

SNR radius 

22 

V s (km/s)* 

Shock velocity 

500 

* (yr)* 

SNR age 

17400 

n 0 (cm 3 ) 

Average ambient density 

0.12 

E 51 (10 51 erg) 

Explosion energy 

0.5 

Me) (M e ) 

SN ejecta mass 

4.5 

PWN Properties: 

Rpwn (pc)* 

PWN radius 

5.0 

Vp (km/s) 

Pulsar velocity 

400 (north) 

L x ( 2 - 10 ) (erg/s) 

PWN X-ray luminosity 

7:2 10 44 

E'o (erg/s) 

Initial spin-down luminosity 

2:8 10 38 

n 

Pulsar braking index 

3.0 

0 (yr) 

Spin-down timescale 

2000 

B ( G) 

PWN magnetic field 

11 

Density Gradient: 

X 

Density contrast of 12.5 

1.08 

H (pc) 
Orientation 

Characteristic length scale 

5.2 

East /West 


Note. — The SNR properties were derived using the Sedov 
model and the best-fit thermal parameters, assuming a distance 
of 9.0 kpc (Section 5). The magnetic field B was estimated from 
the model in Section 7, while the pulsar velocity v p was estimated 
in Temim et al. (2009). Parameters marked with an asterisk are 
not inputs to the HD simulation. The ejecta mass M e j , pulsar 
braking index, initial spin-down luminosity E o , and the spin-down 
timescale o were adjusted to produce the desired PWN morphol- 
ogy and SNR/PWN dimension at the derived SNR age of 17000 
yr. The density gradient function is described by equation 14 of 
Blondin et al. (2001). 

to n = 3.0, and found an initial spin-down luminosity of 
Eo — 2.8 x 10 38 ergs -1 and an initial spin-down timescale 
of the pulsar tq = 2000 yr. This leads to a current 
spin-down luminosity E = 3.1 x 10 36 ergs -1 . The ob- 
served non-thermal luminosity of the PWN in G327. 1-1.1 
is Lx (2 — lOkeV) = 7.2 x 10 34 ergs - 1 , for a distance of 9 
kpc, which gives E/Lx — 0.023. This ratio is within the 
range observed for other pulsars (Possenti et al. 2002). 
The total SN ejecta mass used for the final run is 4.5 Mq. 
The outputs of the simulation are the total pressure and 
density maps at each time step, as well as maps that 
separately trace the ISM, SN ejecta, and PWN material. 
We also trace the age of the particles injected by the pul- 
sar in order to determine where in the PWN we should 
expect spectral steeping due to synchrotron aging. This 
can then be compared with the spectral steepening that 
we measure from the PWN spectra in Section 5. 

6.2. Comparison with G327. 1-1.1 

The results of the HD simulation at four different ages 
are shown in Figure 6. The combined effect of the higher 
ambient density in the west and the northward displace- 
ment of the PWN due to the pulsar’s motion result in 
an asymmetric interaction of the RS with the PWN from 
the NW direction. The interaction begins at an age of 
3000 yr, with the RS reaching the entire surface of the 
PWN at around 7000 yr. At this stage, the Rayleigh- 
Taylor instabilities at the PWN/RS interface also be- 


come more prominent. The RS sweeps the PWN away 
from the pulsar and displaces it eastward of the explo- 
sion center. Once the RS reaches the pulsar itself, at an 
age of ~ 12000 yr, a narrow finger/trail of PWN mate- 
rial begins to form. Since the pulsar continually injects 
wind-producing particles as it moves to the north, in the 
absence of any density gradient, the PWN trail would 
be oriented northward. However, in our case, the north- 
ward motion of the pulsar combined with the preferential 
arrival of the RS from the west causes the trail to be ori- 
ented along the NW direction. 

The simulation suggests that the displacement of the 
relic PWN in G327. 1-1.1 is mostly determined by the 
orientation of the ambient density gradient, while the 
orientation of the trail behind the pulsar is determined 
by the combination of the density gradient and pulsar’s 
velocity. The formation and thickness of the trail also 
strongly depends on the explosion energy and the pul- 
sar’s spin-down luminosity. A spin-down time allowing 
for an initial spin-down luminosity below 10 37 ergs -1 will 
generally struggle to displace the incoming shock and 
form a narrow trail, while a higher energy output will 
easily maintain a larger bubble around the pulsar. None 
of the simulation runs produced non-thermal structures 
ahead of the pulsar that might explain the prongs and 
arcs observed in the Chandra images. 

At an age of 17000 yr, the Sedov estimated age for 
G327. 1-1.1, the morphology closely resembles the radio 
and X-ray observations of the SNR. The physical sizes of 
the SNR and PWN radii correspond well to those mea- 
sured for G327. 1-1.1 assuming a distance of 9 kpc. The 
simulation also shows that the SN ejecta mixes with the 
PWN material during the RS interaction, as suggested 
by the observed morphology of the soft X-ray emission 
in Figure 1. 

The simulation also keeps track of the age of the par- 
ticles injected by the pulsar. Figure 7 shows the map 
of the particle age parameter at the present SNR age of 
17000 yr. The freshly-injected PWN material is shown 
in blue, and the older particles in the light shades of red. 
As expected, the morphology of the aged particles in red 
resembles the radio emission from the PWN in G327.1- 
1.1, while the distribution of younger particles in blue 
and green traces the X-ray emission along the trail. Our 
analysis of the X-ray spectroscopy shows a steepening in 
the X-ray trail of Ar = 0.52 ±0.17 (Section 5), consistent 
with the steepening expected from synchrotron aging of 
particles as they travel away from the pulsar. Using a 
magnetic field of 11 pG (see Section 7), we estimate that 
the synchrotron lifetime over which the steepening is ex- 
pected to occur is T S yn = 1700 yr. This means that in 
order to observe a spectral steeping of AT = 0.5 along 
the trail, the particles would have to had traveled for at 
least 1700 yr. Figure 7 shows that the particles at the 
base of the trail are ~ 4000 yr old, consistent with the 
scenario in which the observed steeping in the X-ray trail 
is caused by synchrotron aging. 

7. BROADBAND EMISSION 

To investigate the broadband emission from G327.1- 
1.1, we have used the semi- analytical model of Gelfand 
et al. (2009) for the evolution of a PWN within the con- 
fines of an SNR. This one dimensional model simulta- 
neously treats the magnetic, dynamical, and radiative 



DEEP CHANDRA OBSERVATION OF G327. 1-1.1 


11 


Age (yr) 



Fig. 7. — Map of the PWN from the HD simulation described in 
Section 6, showing the age of the particles injected by the pulsar 
at an SNR of 17000 years. The region shown is the same as the 
blue colored PWN region of the lower right panel of Figure 6. The 
freshly-injected particles (blue) give rise to the X-ray emission in 
the cometary structure and the trail, while the aged PWN material 
shown in red corresponds to the morphology of G327. 1-1.1 observed 
at radio wavelengths. 

evolution of the system, and given a set of input pa- 
rameters, predicts properties such as the magnetic field, 
PWN and SNR sizes and expansion velocities, and the 
broadband spectrum of the PWN, as a function of the age 
of the system. We use an ambient density of 0.12 cm -2 
based on our X-ray estimates for a distance of 9 kpc, 
and evolve the system until the SNR radius matches the 
observed value of 22 pc. We use the results of our HD 
modeling for the explosion energy, ejecta mass, initial 
pulsar spin-down power, and pulsar spin-down timescale 
(see Table 3). We assume the injection of a magnetized 
electron/positron wind with a broken power law for the 
particle spectrum. The spectral index for the low en- 
ergy component is fixed at 1 .48 based the observed radio 
index (Ma et al. 2015), while that for the high energy 
component is fixed at 2.2 based on our measured X-ray 
spectral index for the emission in the compact source 
immediately surrounding the pulsar (see Section 5). 

The resulting evolution of the SNR and PWN radii 
is shown in Figure 8. The model reaches the radius of 
G327. 1-1.1 (indicated by the blue dashed line) at an age 
of ~ 17800 yr. By this time, the PWN has been com- 
pressed by the SNR reverse shock, and has begun to 
expand again, as seen from the evolution of the PWN 
radius. At the current age, the PWN radius is ~ 5 pc, 
in good agreement with the size of the relic nebula in 
G327.1-1.1. 

The magnetic field in the PWN, shown in the upper 
panel of Figure 8, is evolved through the process of ex- 
pansion and RS compression under the assumption that a 
constant fraction r / b of the spin-down power is injected 
in the form of magnetic flux, with the remaining frac- 
tion rje representing the injected particle flux. Radiative 
and adiabatic losses are calculated throughout the evolu- 
tion, and the final spectrum is obtained from the result- 
ing electron population using the evolved magnetic field 
to calculate the synchrotron radiation. Inverse-Compton 
emission is calculated for scattering of the electron popu- 
lation from photon fields associated with the CMB, local 
starlight, and infrared emission from local dust (Strong 
et al. 2000). We find that an electron spectrum com- 
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Fig. 8. — Top: Model for emission from the evolved PWN in 
G327. 1-1.1 (see Section 8) along with observed radio emission from 
the PWN in black (Ma et al. 2015), X-ray flux and spectral index in 
cyan (Section 5), H.E.S.S. -ray emission in magenta (Acero et al. 
2012) and Fermi upper limits in red (Acero et al. 2013). Bottom: 
Time evolution of the PWN magnetic field (upper panel) and SNR 
and PWN radii (lower panel) for G327. 1-1.1. The observed SNR 
and PWN radii are indicated by the blue and red horizontal dashed 
lines, respectively, and the age at which these radii are reached in 
the model is indicated by the green vertical dashed line. 


prising 99.7% of the spin-down power at injection, with 
an energy break at ~ 300 GeV, produces a good rep- 
resentation of the broadband spectrum, which is shown 
in Figure 8. Here the X-ray spectrum (cyan) represents 
that of the entire PWN, as derived in Section 5, and the 
H.E.S.S. spectrum (magenta) is that reported by Acero 
et al. (2012). The Fermi LAT upper limits (red) are 
taken from Acero et al. (2013). The modeled magnetic 
field at the current SNR age is 10.8 /iG, which based on 
the ~ 300 GeV intrinsic break in the electron spectrum, 
predicts a spectral break in the mid-infrared band. 

We note that the treatment here is oversimplified in 
several ways as a result of the spherically-symmetric 
evolution into a constant density in the semi-analytical 
model. This ignores the effects of the asymmetric crush- 
ing of the PWN, so the resulting magnetic field evolution 
is not likely to reflect the more complicated behavior ex- 
pected in G327. 1-1.1. Despite these inadequacies, the 
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model provides a reasonable representation of the broad- 
band spectrum, and of the PWN radius - at least for the 
relic nebula. While the spectrum of the high energy elec- 
trons corresponds to that derived from the emission close 
to the pulsar, the overall X-ray spectrum is produced by 
electrons near the maximum energy, which are suffering 
from loss-induced steepening of the electron spectrum. 

8. EVOLUTIONARY STATE 

The modeling described in Sections 6 and 7 indicates 
that the morphology and broadband spectrum of G327.1- 
1.1 is well described by the evolution of a composite SNR 
whose PWN has been crushed by the RS. A scenario in 
which the SN blast wave expands in non-uniform ISM 
density that decreases from west to east, and in which 
the pulsar is moving to the north with a velocity of 400 
kms -1 , matches the observed morphology of the G327.1- 
1.1 at radio and X-ray wavelength remarkably well. The 
observed structure of the SNR shell and relic PWN at 
radio wavelengths, as well as the morphology of the X- 
ray PWN and trail, are reproduced by the HD simula- 
tion at the age of ~ 17000 yr, consistent with the age 
of G327-1.1 as derived from the Sedov model. Despite 
the oversimplification of spherical symmetry, the semi- 
analytic model of the evolution of G327. 1-1.1 can match 
the broadband spectrum of the PWN emission from radio 
to 7-ray wavelengths, with input parameters consistent 
with the observed parameters and those derived by the 
HD simulation. 

In earlier work, we speculated that the cometary X- 
ray PWN in G327. 1-1.1 is likely not a bow shock, due 
to the fact that we observe emission ahead of the pul- 
sar (Temim et al. 2009). In Section 4.1.1, we show that 
the compact source immediately surrounding the pulsar 
shows evidence for a torus-like structure embedded in 
the cometary PWN. This provides further evidence that 
the newly-formed PWN is not in the bow-shock stage 
of evolution, but that its cometary morphology is most 
likely a result of the passage of the RS that reaches the 
PWN from the NW direction (see Section 6) . In addition, 
MHD simulations of pulsar bow-shock and subsonic neb- 
ulae show that in a fully subsonic model, the cometary 
nebula and X-ray trail are more confined and have typ- 
ical particles ages on the order of a few hundred years 
(Bucciantini et al. 2005; Swartz et al. 2014). This is 
consistent with our observations of the size and spectral 
steepening of the X-ray PWN in G327. 1-1.1. We note 
that these same models predict a much smaller size for 
the termination shock radius than the size of the com- 
pact source that we attribute to a torus-like structure, 
suggesting that the structure that we are observing is a 
Crab-like torus that persists downstream of the termi- 
nation shock, rather than a ring from the termination 
shock itself. 

While the above scenario provides striking consistency 
with the observed spectral and spatial properties of 
G327. 1-1.1, some outstanding questions still remain. The 
prong and arc structures are the most puzzling and diffi- 
cult to explain, especially since they don’t seem to orig- 
inate at the pulsar itself. None of the HD simulations of 
the PWN crushing process produced non-thermal emis- 
sion ahead of the pulsar’s motion. One possibility is that 
these are structures associated with enhancements in the 


magnetic field. In Section 5, we noted that the prong and 
arc structures are reminiscent of a funnel-shaped outflow 
aligned with the X-ray trail (i.e., the radio finger). The 
fact that this “outflow” is oriented in the NW direction, 
the effective direction from which the RS first reached 
the pulsar and the PWN, may suggest an association 
with the RS/PWN interaction. Three-dimensional HD 
simulations of this system may provide additional clues 
into the origin of these unusual structures. 

9. CONCLUSIONS 

In this paper, we analyzed 350 ks Chandra observa- 
tions of the composite SNR G327.1-1 and showed that 
the emission originates from a thermal shell with a tem- 
perature of 0. 29^0:04 keV, and a PWN with an aver- 
age photon index of 2.15 ^q: 04- The Chandra imaging 
reveals the complex morphology of the PWN; a compact 
source surrounding the presumed pulsar, embedded in a 
cometary X-ray nebula with a trail of emission behind it. 
The deeper X-ray observation reveals non-thermal emis- 
sion coincident with the radio PWN, with a steep pho- 
ton index of 2.581q!io- There is also evidence for thermal 
emission with enhanced abundances in the PWN region. 
Morphological evidence suggests that it originates from 
SN ejecta mixed-in with the PWN material following 
a RS interaction. Unusual prong-like structures origi- 
nate from the vicinity of the PWN and extend into large 
arcs. The structure is reminiscent of a funnel-shaped 
outflow towards the NW direction. The deconvolution of 
the compact source immediately surrounding the pulsar 
shows some evidence for a torus-like shape, which sup- 
ports our conclusion from previous work that the PWN 
has not formed a bow shock and is expanding subsoni- 
cally. The average absorbing column density for the SNR 
is approximately 2 x 10 22 cm -2 , with evidence of higher 
values in the western portion of the SNR. 

HD simulations of a composite SNR expanding in a 
non-uniform density medium and hosting a rapidly mov- 
ing pulsar reproduces the morphology of G327. 1-1.1 ex- 
ceptionally well. The simulation was constrained by the 
observationally determined SNR and PWN properties 
and evolved to an age of ~ 17000 yr, the estimated age 
of G327. 1-1.1 from the Sedov calculation. It shows that 
the morphology of G327. 1-1.1 most strongly depends on 
the strength and orientation of the ISM density gradi- 
ent, and the pulsar’s velocity, spin-down luminosity, and 
spin-down time frame. A scenario in which the ambi- 
ent density gradient decreases from west to east, and the 
pulsar moves to the north with a velocity 400 kms -1 , 
causes the RS to crush the PWN preferentially from the 
NW direction and produces a morphology that closely 
resembles the radio and X-ray observations of the SNR. 
The RS/PWN interaction scenario can also reproduce 
the broadband spectrum of the PWN at radio, X-ray, 
and 7-ray bands. Future modeling of similar systems 
whose physical properties can be constrained by observa- 
tions promise to significantly improve our understanding 
of the evolution of composite SNRs. 
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